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Abstract  -  In  this  paper  we  report  our  experiences  with  a  new 
texture  coarseness  measure,  as  a  step  towards  automation  of 
metaphase  finding  in  cell  proliferation  studies.  This  measure  is 
highly  specific  to  grey-level  inter-chromosome  coarseness 
features  in  microscopic  images  of  metaphase  spreads,  and  allows 
to  quantify  the  texture  of  the  cytological  objects  analysing  the 
intensity  profile  between  chromosome-extrema  samples. 
Chromosome  fragments  produce  patterns  of  pixels  at  low- 
resolution,  and  the  local  neighbourhood  of  their  individual 
extrema  presents  a  characteristic  coarseness,  along  intensity 
profiles  on  random-oriented  test  lines.  Results  of  its  use  on 
images  of  fields  of  metaphases  and  artefacts  are  compared  with 
some  representative  texture  measures,  outperforming 
metaphase  detection  and  artefacts  elimination.  This  coarseness 
feature  provides  a  specific  metaphase  signature  that  can  be  used 
in  conjunction  with  other  morphological  and  textural 
parameters  for  automated  metaphase  discrimination. 
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extrema. 

I.  Introduction 

The  microscope  discrimination  of  metaphase  spreads  is  a 
method  widely  and  routinely  used  for  the  detection  of 
chromosomal  aberrations  in  medical  genetics,  toxicological 
studies  of  mutagenicity,  cell  proliferation  evaluation  for 
drugs  safety,  and  radiation  health  monitoring,  amongst 
others. 

The  increasing  demands  on  health  services  to  provide 
genetic  counselling,  pregnancy  screening,  evaluation  of  drugs 
in  cellular  proliferation  and  radiation  health  monitoring  have 
added  to  the  pressures  on  cytogenetic  laboratory  services. 
The  task  involved  in  locating  a  number  of  individual 
metaphases  is  very  time-consuming,  even  at  low  resolutions. 
In  both,  routine  and  research  environments,  the  involved  task 
are  tedious  and  tiring,  but  the  importance  and  consequences 
of  the  results  of  screening  process  demand  consistency  and 
accuracy  from  highly  trained  and  qualified  staff.  The 
combined  effect  of  these  factors,  coupled  with  the 
omnipresent  need  to  avoid  further  cost,  has  prompted  and 
accelerated  the  design  and  development  of  systems  for 
automated  metaphase  spread  location. 

Several  approaches  and  systems  known  as  metaphase 
finders  have  been  proposed.  The  first  instrument  for  semi- 
automated  metaphase  finding  was  described  in  a  paper 
published  in  1963  [1].  This  system  was  initially  concerned  to 
locate  metaphase  and  was  based  on  a  conventional  optical 
microscope  with  a  scanning  stage  using  video  imaging.  Other 
attempts  to  semi-automate  metaphase  finding  using  a  closed- 
circuit  TV  camera  were  reported  in  [2,3].  More  recent  works 


report  the  use  of  commercial  metaphase  finders  and  their 
performance  [4,5,6].  Although  significant  progress  has  been 
made  in  automatic  metaphase  finding  technology,  the  total 
time  spent  in  the  analysis  of  slides  and  a  greater  accuracy  still 
represent  important  challenges  to  biomedical  engineers. 

In  the  present  study,  we  concentrate  our  efforts  in 
extracting  a  texture  measure  from  the  inter-chromosome 
profile  information,  achieving  results  that  substantially 
enhance  metaphase  recognition  and  artefact  discrimination. 
The  texture  measure  proposed  was  tested  in  a  current  cell 
proliferation  study.  Working  in  a  tenfold  magnification  and  a 
visual  field  of  610x477pm,  this  profile  feature  allows  the 
identification  of  variable-shaped  metaphases  in  the  presence 
of  many  artefacts.  The  identification  performance  was 
compared  with  that  of  classical  texture  features  such  as  RED 
[7]  and  those  derived  from  grey-level  co-ocurrence  matrices: 
contrast,  entropy,  energy,  angular  second  moment  [8] and 
MDWRE  [  9]. 

II.  Metaphase  textural  information 

The  human  expert  identifies  metaphase  spreads  mainly  by 
observing  the  internal  texture  produced  by  the  chromosomes 
inside.  In  this  sense,,  the  microscope  images  show  that  texture 
in  metaphase  is  composed  of  cytological  material  and  image 
noise,  but  it  also  comprises  digitalized,  sub-sampled  and 
overlapped  chromosome  fragments.  This  composition 
organizes  grey-level  variations  and  specific  pits 
(chromosome  fragments)  into  characteristic  signatures  along 
intensity  profiles  on  random-oriented  test  lines.  Our 
experience,  and  that  of  other  researchers,  has  shown  that 
relationships  between  extrema  features  are  meaningful 
characteristics  of  such  signatures  [10].  Moreover,  the  semi- 
structural  character  of  such  textures  (with  the  sub-sampled 
chromosomes  constituting  the  texels )  makes  them  suitable  for 
run-length  analysis  [11],  and  stereological-feature  extraction 
by  profile  analysis  along  test  lines  [12]. 

Grey-level  images  I(x,y)  are  often  studied  as  surfaces,  with 
intensity  seen  as  height,  and  texture  appearing  as  relief  whose 
extrema  features  combine  intensity  and  shape  information. 
We  propose  to  simplify  2D  feature  analysis  by  interpreting 
the  statistical  features  of  ID  profile  samples  between  extrema 
dots.  The  mean  and  the  standard  deviation  of  intensity  height 
along  such  profiles,  that  is,  the  linear  RMS-roughness  value, 
cannot  serve  alone  to  discriminate  different  signatures.  In  the 
other  hand,  a  2D  RMS-roughness  evaluation  averages-out  all 
features,  while  in  the  present  study  we  seek  inter-extrema 
texture  signatures,  which  are  best  evaluated  by  analysing  just 
ID  profiles  between  extrema  points,  and  considering  the 
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ranges ,  that  is,  the  largest  differences  between  maximal  and 
minima. 


III.  Inter-extrema  profile  texture  feature  (iepr) 


In  a  profile  or  region  containing  N  height  (grey  level) 
samples,  their  values  can  be  ranked  to  obtain  the  n-order 
statistics ,  with  n=l,...,  N.  Besides  traditional  RMS-roughness 
measures,  statistical  rank  filters  like  the  median,  or  the 
trimmed  mean  [10]  can  be  used  to  obtain  simple  coarseness 
descriptors.  In  particular,  the  median  rank  filter  is  the  N/2- 
order  statistics,  if  N  is  even,  and  average  of  the  (N-l)/2  and 
(N+l)/2  values  for  N  odd,  and  it  constitutes  a  robust  RMS- 
roughness  estimator.  These  observations  prompted  us  to 
focus  on  the  very  distributions  of  ranges  (the  largest 
difference  between  maximal  and  minima  in  each  profile),  and 
device  a  simple  coarseness  feature.  While  the  median  rank 
filter  (a  “robust”  RMS  roughness  estimator)  is  representative 
of  height  differences,  coarseness  is  better  described  by  the 
distribution  of  larger  height  differences.  Thus,  the  histogram 
median  (not  the  rank  filter  median)  is  used  once  all  profile 
ranges  (the  N-order  statistics)  are  counted.  Besides 
extracting  the  median  of  the  range  population,  we  select  only 
those  profiles  (1)  between  chromosome  intensity  minima  (at 
each  chromosome  pattern),  and  (2)  within  a  window 
representing  the  inter-extrema  mean  separation.  Note  that 
such  profiles  are  equivalent  to  those  of  diameter  samples 
across  the  neighbourhood  of  each  minimum.  Our  tests  in  fact 
show  no  difference,  provided  that  the  neighbourhood  size  is 
the  same  as  the  window  size  representing  the  inter-extrema 
mean  separation.  We  formalize  in  mathematical  terms  the 
above  description  as  follows. 

Let  fp  be  the  set  of  N  points  Pt  (xifyt),  with  z=l,...,N,  and 
whose  integer  coordinates,  xhyh  lie  within  the  analysed  object 
boundary.  These  points  PL  are  actually  the  coordinates  of  the 
extrema  (chromosome  or  artefacts  dots)  of  an  intensity- 
surface  image  I(x,y),  see  Fig.  lb.  We  choose  pairs  (Pt ,  Pj  ), 
zV j  at  random,  within  the  window  intervals. 

To  define  the  profile  samples  (Fig.  lc)  along  test  lines,  let 
Pij  ~{^t  =  (1 -  0 P[  +  tPj ,  t  e  [0,1] }  be  the  straight  line 

between  the  points  Pt  and  Pj  ,  and  let  I{.  (£  t )  be  the  grey 

level  at  £t.  then: 


\l{Pi)  if  t  =  0 

\l(Pj)  if  t  =  1 


(1) 


are  the  endpoint  intensities.  Note  that  in  a  discrete  space,  I(xt, 
yt)  may  not  be  defined.  In  such  a  case,  its  value  is  estimated 
by  a  bilinear  interpolation,  from  the  nearest  neighbours. 

Let 


be  the  range  of  the  profile-intensity  defined  by  the  straight- 
line  £u  where  tE  [0,1].  The  coordinates  of  the  initial  point  Pt 
and  the  final  point  Pj  of  the  profile  are  within  a  rectangular 

window  where  \xt  -  Xj  |  <  77  and  | yt  -  y j  |  <  17 ;  in  our  case 

77=3.  Since  chromosome  intensities  vary  randomly,  their 
extrema  location  and  intensities  are  random  variables,  and  so 
their  difference. 

Let  H(  dif  ij  )  be  the  histogram  of  all  ranges,  defined  in  the 
equation  (  2  ),  over  the  set  p,  with  grey  level  bins  in  [0,255]. 
Let  H  be  the  median  of  the  histogram  //,  then  we  define 

IEPR  =  ///(/max  -  min{  min  Iy(£t)})  ,  (  3  ) 

/  hj 

where  Imax  is  the  maximum  saturation  grey  level  of  the  image 
and  £e[0,1].  Thus,  IEPR  is  a  texture  measure  of  local 
coarseness,  along  the  intensity  profile  between  extrema. 
Normalization  allows  proper  comparison  under  different 
contrast  condition. 

Denominator  7max  -min {min  /■■  (£ t)}  is  never  null,  since 

ij  itep9 

profile  endpoints  include  always  the  chromosome  (or 
artefact)  minima. 


Distance  (Pixel) 


Fig.  1.  Inter-chromosome  (or  inter-extrema)  profile  analysis,  a) 
metaphase  spread,  b)  detected  extrema  in  a  horizontal  and  vertical  scanning 
of  figure  1  (a),  c)  profile  intensity  graphic  between  two  extrema  points  Pi  and 
pj- 


IV.  Experimental  results 


difij  =  max  Iy  {lt)~  min  Iy (£ t ) 


To  show  the  performance  of  IEPR,  as  a  discriminating 
(  2  )  feature  in  the  classification  of  metaphases,  we  have  applied  it 
to  the  problem  of  metaphase  finding  in  cell  proliferation 
studies  [13]. 


A.  Test  sets 

Twenty-two  slides  from  different  experiments,  and  six 
donors  where  randomly  chosen  to  create  two  test  sets, 
composed  of  76  metaphases  and  191  artefacts. 

Each  test  image  was  digitised  on  a  512x480-pixel  grid  with 
grey-level  resolution  of  256  levels. 

It  is  important  to  remark  that  artefacts  on  slide  image 
present  wide  variation  in  shape  size  and  texture.  This 
indeterminacy  makes  artefacts  to  have  a  high  possibility  of 
being  a  source  of  errors  (false  positives,  mainly)  in  the 
automatic  metaphase-fmding  problem.  Moreover,  the 
probability  of  artefacts  is  extremely  high  (one  metaphase  for 
a  thousand  of  artefacts  in  mitotic  index  analysis).  Fortunately, 
a  large  percentage  of  these  artefacts  are  outside  the  size  range 
of  metaphases  we  want  to  score,  and  it  is  then  possible  to 
apply  an  outlier  exclusion  criterion.  Nevertheless,  the 
remaining  artefacts  are  the  most  difficult  to  discriminate,  as 
many  of  them  may  be  similar  to  the  objects  of  interest.  The 
metaphase  correctly  recognized  represents  the  true  positive 
metaphase  subset,  denoted  as  TP. 

B.  Evaluation  of  IEPR  for  the  test  sets 

IEPR  was  calculated  for  all  objects  from  the  two  test  sets 
(metaphase  spreads  and  artefacts).  Results  are  summarized  in 
Table  1.  The  nuclei  and  artefacts  objects  miss-classified  as 
metaphase  represents  the  false-positives  metaphase  subset, 
denoted  by  FP.  Comparison  of  performances  included  IEPR, 
MDWRE,  relative  extrema  density  (RED) [7]  and  the  most 
meaningful  measures  from  grey-level  co-occurrence  matrices 
(0°  angular  relationship):  contrast,  energy,  entropy  and 
angular  second  moment  (ASM)  [8]. 

Table  i. 

IEPR  PERFORMANCE  COMPARISON  WITH  OTHER  TEXTURE  FEATURES. 


PERFORMANCE=(TP-FP)/  TOTOBJ 


TP 

metaphase 

FP 

metaphase 

Performance 

IEPR 

96% 

5% 

0.9079 

MDWRE 

84% 

14% 

0.7105 

Contrast 

96% 

25% 

0.6447 

RED 

74% 

29% 

0.4473 

Energy 

32% 

48% 

0.0263 

ASM 

33% 

65% 

-  0.2895 

Entropy 

18% 

75% 

-  0.3684 

V.  Discussion  and  conclusions 


We  have  proposed  and  evaluated  in  this  work  an  inter¬ 
chromosome  specific  measure  on  profile  range  distributions, 
used  as  a  characteristic  signature  of  metaphase  spreads,  as 
compared  with  present  artefacts.  This  signature  relates  local 
range  distributions  with  the  local  coarseness  characteristically 
produced  by  chromosome  fragments. 

Based  on  the  percentage  of  true  positives  (TP)  and  false 
positives  (FP),  a  comparison  between  IEPR  and  other  texture 
feature’s  performances  is  achieved.  The  performance  of 
IEPR,  as  a  texture  feature,  surpasses  by  22%  the  best  of  the 


compared  features  (MDWRE)  for  metaphase  finding  (Table 

1). 

Variograms,  run-lengths  and  co-occurrence  features 
describe  better  spatial  dependency  and  extrema  separation 
relationships,  while  intensity  signatures  describe  the  kind  of 
coarseness  differences  observed  among  stimulated  nuclei, 
artefacts,  cluttered  or  dispersed  chromosomes.  Other 
approaches  (e.g.  multi-resolution,  Markov  random  fields  and 
Gabor  filters)  require  higher  resolution  and  precious 
processing  time  to  give  meaningful  discrimination  results. 

The  proposed  coarseness  measure  represents  a  reliable  and 
cost  effective  metaphase  signature  that  can  be  used  in 
conjunction  with  other  morphological  and  textural  parameters 
for  automated  metaphase  discrimination. 

Acknowledgements 

This  work  was  supported  by  the  Programa  de  Apoyo  a 
Proyectos  de  Investigacion  e  Innovacion  Tecnologica 
(PAPIIT-  INI 05 500)  to  GCB.  We  thanks  Blanca  Taboada  her 
valuable  technical  support. 

References 

[1]  Ledley  R.,  Ruddel  F.  and  Migeon  B.  (1963),  Proc.  16th  Ann. 
Conf.  Eng.  Med.  Biol. 

[2]  Johnson  E.  and  Goforth  L.  (1974).  “Metaphase  spread  detection 
and  focus  using  closed  circuit  television”,  J.  Histochem.  Cytochem., 
vol.  22, No.  7,  pp.  536-545. 

[3]  Green  D.  and  Neurath  P.  (1974).  “The  desing,  operation  and 
evaluation  of  a  high  speed  automatic  metaphase  finding”.  J. 
Histochem.  Cytochem.,  vol.  22, No.  7,  pp.  531-535. 

[4]  Castleman  K.R.(1992),  “The  PSI  automatic  metaphase  finder”,  J. 
Radiat.  Res.,  33,  (Supl.),  pp.  124-128. 

[5]  Weber  J.,  Scheid  W.  and  Traut  H.  (1992).  “Time-saving  in 
biological  dosimetry  by  using  the  automatic  metaphase  finder 
Metafer2”,  Mutat.  Res.  Vol  272,  No.  1,  pp-31-34. 

[6]  Piper  J.,  Poggensee  M.,  Hill  W.,  Jensen  R.,  Ji  L.,  Poole  I.,  Stark 
M.  and  Sudar  D.  (1994).  “Automatic  fluorescence  metaphase  finder 
speeds  translocation  scoring  in  FISH  painted  chromosomes”, 
Cytometry  ,  Vol.  16,  No.l,  pp.  7-16 

[7]  Rosenfeld  A.  and  Troy  E.  (1970).  “Visual  Texture  Analysis”, 
Technical  Report,  University  of  Maryland,  College  Park,  pp.70-1 16. 

[8]  Haralick  R.  (1979),  “Statistical  and  structural  approaches  to 
texture”,  Proc.  IEEE,  vol  67,  pp.  786-804. 

[9]  Corkidi  G.,  Vega  L.,  Marquez  J.,  Rojas  E.,  Ostrosky-Wegman  P. 
(1998),  “Roughness  Feature  of  Metaphase  Chromosome  Spreads 
and  Nuclei  for  Automated  Cell  Proliferation  Analysis”,  Med.  Biol. 
Eng.  Comput.,  vol  36,  pp.  679-685. 

[10]  Haralick  R.,  and  Shaphiro  L.  (1992)  “Computer  and  Robot 
Vision”,  Volume  1.  Addison  Wesley. 

[11]  Galloway  M.  (1975),  “Texture  Analysis  Using  Gray  Level  Run 
Lengths”.  Computer  Graphics  and  Image  Processing  No.  4,  172- 
179. 

[12]  Marquez  J.  (1999),  “Analysis  and  Visualization  of  Three 
Dimensional  Complex  Structures:  a  Boundary  Based  Approach”. 
PhD  Thesis,  Dept.  TSI,  ENST,  pp.  197-210. 

[13]  Rojas  E.,  Montero  R.,  Herrera  L.A.,  Sordo  M.,  Gonsebatt  M.E., 
Rodriguez  R.,  Ostrosky-Wegman  P.  (1992),  “Are  mitotic  index  and 
lymphocyte  kinetic  reproducible  endpoints  in  genetic  toxicology 
testing?”,  Mutation  Res.,  282,  pp. 283-286 


